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The variant of the single-impurity Kondo problem in which the conduction-band density of states has 
a power-law pseudogap at the Fermi energy is known to exhibit a zero-temperature phase transition 
at a finite exchange coupling. The critical properties of this transition are studied both for N — 2 
and for N S> 1, where N is the spin degeneracy. The critical exponents are consistent with a 
simple scaling form for the free energy. For any finite N, the temperature exponent of the local spin 
susceptibility at the critical Kondo coupling varies continuously with the power of the pseudogap. 
This raises the possibility that a single-particle pseudogap is responsible for the anomalous behavior 
of certain heavy-fermion metals close to a magnetic quantum phase transition. 

PACS numbers: 71.10.Hf, 75.20.Hr, 71.27. +a, 71.28.+d 



A number of heavy-fermion metals exhibit non- Fermi- 
liquid behavior which seems to arise from proximity to a 
zero-temperature magnetic ordering transition jjj . These 
systems represent a puzzle in that temperature exponents 
describing various transport and thermodynamic proper- 
ties differ from the predictions of the standard Ginzburg- 
Landau-type theory for spin fluctuations ||||]. Recent 
neutron scattering experiments in CeC^-rcAu^ suggest 
that this discrepancy may arise from an unusual mo- 
mentum and energy dependence of the spin susceptibility 
In particular, an anomalous energy dependence ex- 
tends over a wide range of momenta, suggesting that the 
local spin susceptibility is also anomalous [j5j . The micro- 
scopic origin of this local anomaly is unclear at present. 

This paper explores the possibility that an anomalous 
local susceptibility arises from a pseudogap in the single- 
particle spectrum. (In quantum-critical systems which 
display superconductivity, such as CePd2Si2 and Celn3 
||, superconducting fluctuations may generate a pseu- 
dogap j7|. Alternative mechanisms Q] may be at work 
in other systems.) To this end, we study a variant ||] 
of the SU(A)-symmetric Kondo model in which a single 
impurity moment (of spin degeneracy A) couples to a 
conduction band described by the (oversimplified) den- 
sity of states l| 



P(e) = 



Po|e| 




for |e| < 1, 
for lei > 1. 



(1) 



In a metal (r = 0), any antiferromagnetic Kondo coupling 
causes the impurity moment to become quenched at suf- 
ficiently low temperatures ju^. With a power-law pseu- 
dogap (r > 0), by contrast, quenching takes place only 
if the Kondo coupling exceeds a nonzero critical value 
H . This threshold is associated with a zero-temperature 
phase transition, which is the focus of this work. 

Our analysis is based on the observation that the crit- 
ical behavior of the power-law Kondo model reveals it- 
self not in the conventional magnetic susceptibility, but 



rather in the response to a local magnetic field coupled 
solely to the impurity spin. We propose a simple scal- 
ing form for the critical component of the free-energy, 
which in turn leads to scaling relations among the criti- 
cal exponents. Within the large- N approach pl| , p^| , we 
calculate the critical exponents on the strong-coupling 
side, and then apply the scaling ansatz to obtain the re- 
maining exponents. For N — 2, all the exponents are 
computed directly using the numerical renormalization- 
group (NRG) method p3|-p^[, and are found to obey the 
scaling relations to high accuracy. Most importantly, the 
pseudogap fundamentally alters the local spin suscepti- 
bility at the critical Kondo coupling, such that (for A 
finite) the critical exponent varies continuously with r. 
We begin with the Hamiltonian 



k.cr 



H K = } , £kC k(T c k(T 



^0 Q t - 



(2) 



where Jo > is the Kondo coupling, cj^ creates a spin- 
a electron (rx = 1,2, ...,JV) at the impurity site, and 
rl a , (i = 1, 2, . . . , - 1) is a generator of SU(iV). We 
employ a pseudofermion representation of the impurity 
spin operator, S — A -1 a , f^Taa'/a', supplemented 
by a constraint flfa = Q- This paper treats systems 
possessing exact particle-hole symmetry, for which Q = 
A/2; the conduction-band dispersion ek is assumed to 
give rise to the density of states in Eq. ([!]) with < r < 1 . 
More general cases will be described elsewhere. 

To study the magnetic response, we add to Eq. (j^) the 
magnetic field terms 



"mag — / 



(3) 



Here m a = sgn(2rr — N — 1) and H represents the uniform 
magnetic field which couples with equal g factors to the 
impurity and the conduction-electron spins, while h is a 
local field coupling only to the impurity spin flq] . 
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Large-N formulation: The large- N limit of the model 
can be analyzed within slave-boson mean-held theory 
p],p|JT^] . In the path- integral representation of the par- 
tition function, a Lagrange multiplier e/ enforces the 
constraint on the impurity occupancy, and Hubbard- 
Stratonovich held 4> decouples the Kondo interaction 
term. The impurity part of the free energy, obtained 
by integrating out the conduction and pseudofcrmion de- 
grees of freedom, takes the form 



F imp (T,H, h) 



H 



N 

-H,h-> -h 



de n F (e+H/2)S (e,h) 



N 

V 



-Qe f , 



(4) 



where n F is the Fermi function and do is the conduction- 
band phase shift: 



(5o(e, h) = tan 



2 ImG (e + iO+) 



e - e f - h/2 - 4> 2 ReG {e + i0+) 



(5) 



The noninteracting local conduction electron Green's 
function satisfies ImGo(e + i0 + ) = —np(e), and for 
|e| « 1 10, 



ReG (e + iO + ) 



-7rp(e)tan(^) sgne, 
-2p(e)ln|e- 1 |, 



< r < 1, 
r=l. 



The mean-field solution is obtained by minimizing -Fi mp 
with respect to e/ and 4> 2 at zero magnetic field. Particle- 
hole symmetry ensures that ey = 0, while 4> 2 is deter- 
mined by the saddle-point equation 



1 

Jo 



de 



n F {e) (e - £/) ImGp(e) 



2 ReG (e)P 



2 ImGo(e)F 



(G) 



which is exact to 0(1/N). 

The critical coupling J c , extracted by the condition 
that 4> vanishes as Jo — » J+ from Eq. (0) with T = H = 
h = 0, satisfies pqJ c = r ]§). Just above the critical 
coupling, we find that [n7| 



2 oc (Jo - J c ) 1,f -\ f = min(ri). 



Local vs impurity susceptibility: An important clue 
to the nature of the phase transition for r > comes 
from the contrast between the zero-temperature limits of 
the impurity susceptibility Ximp = -d 2 F imp /dH 2 \ H=h=0 
and the local susceptibility xioc = —d 2 F imp /dh 2 \ H=h=0 . 



NRG results [Jl4|,g_5j| for degeneracy N = 2 indicate that 
whereas limT^o^Ximp is nonzero for all values of Jo, 
limj^o Txioc goes continuously to zero as the critical 
coupling is approached from below and liniT^o Txioc = 
for all Jo > J c . We now show that in most cases, the same 
distinction can be made in the large- N limit . 



Starting from Eq. (Q), the impurity susceptibility can 
be written 



N 



Ximp 



1 ^ dn F (e) dS o {e,0) 



47r de de 



FromEq. (§), <5 (e->0,0) w (1 + r sgne)7r/2 @, so 



lim Tximp = Nr/ 16. 



(8) 



(9) 



This large- N result reduces in the case of iV = 2 to the 
value r/8 obtained within the NRG approach pl| . 

The local spin susceptibility can be evaluated directly 
by differentiating Eq. (|4|). Two very different types of 
behavior are predicted, depending on the value of r. For 
< r < i, we find that 



: (J > J c ,T = 0) oc (Jo- J c )- 1/r . 



(10) 



This is consistent with the result obtained for N = 2, 
namely limT^o Jxioc = 0. Note that unlike Ximp, which 
is independent of Jo throughout the strong-coupling 
regime, xioc(T = 0, J ) diverges as J — > J+. 

For i < r < 1, the large- TV result is xioc = oo- This 
unusual behavior has a counterpart in the case N = 2, 
where J c = oo, i.e., the strong-coupling regime is com- 
pletely inaccessible Away from strict particle-hole 
symmetry, a critical xioc is recovered for all r > 0; how- 
ever, the treatment of particle-hole asymmetry lies be- 
yond the scope of the present paper, so we focus hence- 
forth on the range < r < \. 

General scaling analysis: In the standard (r = 0) 
Kondo problem, the zero-temperature phase transition 
between weak and strong coupling has the Kosterlitz- 
Thouless form. Based on the Coulomb-gas picture of the 
partition function [|l0|, one expects that the Kosterlitz- 
Thouless nature of the transition is special to r = 0. For 
r > 0, by contrast, it is natural to assume that the tran- 
sition is continuous (second-order). 

At a continuous transition, the singular component of 
the free energy F s should exhibit scaling. The deviation 
of the Kondo exchange from its critical value, Jo — J c , 
is clearly a relevant variable. The different behaviors of 
Ximp and xioc reported above suggest that the local mag- 
netic field h (rather than the uniform magnetic field H) 
should act as a second scaling variable. F s is therefore 
expected to have the form 



F s =Tf(\J -J c \/T a , \h\/T b ) 



(11) 



Alternatively, one can characterize the critical behav- 
ior by a set of exponents (3, 7, 6, and x, defined as follows: 



AM Jo < Jc, T = 0, h = 0) cx ( J c - J f, 
Xioc( Jo > Jc, T — 0) oc (J - J c )~ 7 , 
Af loc (J = J c ,T = 0) oc H 1 / 5 , 

Xloc(Jo = Jc) oc T~ x . 



(12) 
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The local magnetization M\ oc (h = 0) = —dF s /dh\h=o 
acts as the order parameter for the transition and 
Xioc = —d 2 F s /dh 2 \h = o is the order-parameter suscepti- 
bility. The critical exponents can readily be written in 
terms of a and 6 introduced in Eq. (11): (3 = (1 — b)/a, 
7 = (26 - I) /a, S = 6/(1 - 6), and x = 2b - 1. These 
expressions lead to a pair of scaling relations among the 
critical exponents, e.g., 



6 



(l + x)/(l-x) 
1 {l-x)/{2x). 



(13a) 
(13b) 



Critical behavior in the large-N limit: Equation ( |To| ) 
implies that 7 = 1/r for < r < \. Since the large- 
iV mean- field solution is valid only for Jo > J c , it is 
not possible to determine the exponents x, <5, and (3 di- 
rectly. These exponents may be inferred, however, from 
the general scaling analysis outlined above. Combining 
Eqs. (|M@), we find that 

F S (J > J C ,T = H = h = 0) = -NA(J - J c ) x l\ (14) 

where A is an r-dependent coefficient and f = min(r, i) 
as before. If we assume the validity of Eq. ([TT|), we can 
deduce that a = r and 6 
that the full set of exponents is 



1 for < r < i, and hence 



7= 1/r, 13=1/5 = 0, x = l. 



(15) 



Only 7 shows any r-dependence at the mean-field level. 

We note that —F S (T = 0)/N = Tx, the Kondo scale 
(the inverse correlation length in the temporal direc- 
tion). The vanishing of Tk in power-law fashion as Jo 
approaches J c [see Eq. (jl^)] is precisely the behavior ex- 
pected at a second-order phase transition, and confirms 
that the critical point for r > is not of the Kosterlitz- 
Thouless type. 

Critical behavior for N = 2: We treat the case N = 2 
using a generalization of Wilson's NRG method to de- 
scribe a nonconstant density of states [j^. The local 
magnetization and local susceptibility can be computed 
as M loc = (S z ) and xioc = lhn h -^o(S z )/h, where (S z ) is 
the expectation value of the z component of the impurity 
spin. We have performed calculations for several different 
values of r between 0.1 and 0.45. (As mentioned above, 
the strong-coupling regime disappears for r > |.) In each 
case Mice and xioc behave as described in Eqs. (|T2]), es- 
tablishing the continuous nature of the phase transition. 

Table | lists the exponents extracted from the NRG 
runs along with their estimated random errors. Data at 
the critical point exhibit asymptotic power-laws over at 
least five decades of T or h, resulting in rather precise 
determinations of x or 5, respectively. The determina- 
tion of the other two exponents is hindered by numerical 
rounding errors which eventually cut off the power-law 
variation of M\ oc and xioc on approach to the critical 
coupling. Some of the (3 and 7 values in Table H have 
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P 
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1/5 


1-x 


0.1 


9 






0.005651(3) 


0.0112(1) 


0.15 
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0.095(5) 


7.14(9) 


0.01367(3) 


0.0270(1) 


0.2 


9 


0.153(2) 


5.81(5) 


0.02636(4) 


0.0515(1) 


0.3 
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0.35(1) 


4.39(2) 


07405(51 


1 378(51 
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0.13715(5) 


0.4 


9 


0.90(1) 


3.90(2) 


0.1850(2) 


0.3124(9) 
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3.98(5) 




0.3113(2) 


0.45 


9 


1.91(2) 




0.3135(3) 


0.4771(9) 



TABLE I. Critical exponents for the case N = 2, obtained 
from NRG calculations. Parentheses indicate the estimated 
random error in the last digit. Here A parameterizes the dis- 
cretization of the conduction band. 



been computed from data taken over only one decade of 
\Jq ~~ Jc\- This problem becomes even more pronounced 
as r approaches h, and we have found it impossible to 
make a reliable estimate of 7 for r > 0.4. Very small val- 
ues of r also prove to be inaccessible due to the buildup 
of rounding error during the approach to the limit T = 0. 

Most runs were performed for an NRG discretization 
parameter A = 9, retaining all many-body eigenstates up 
to an energy 50T above the ground state. In an attempt 
to estimate the systematic discretization errors, selected 
runs were performed using A = 3, a value which lies 
closer to the continuum limit (A = 1) but which requires 
much more computer time. It can be seen from Table |l| 
that those exponents computed for both A = 3 and A = 
9 narrowly fail to agree within their estimated random 
errors. We therefore believe that the A = 9 exponents 
provide a fair approximation to the continuum values. 

All four exponents listed in Table | show significant 
r-dependence. In all cases, they sa tisfy the general scal- 
ing relations very well. Equation ( 13a ) relating S and x 
is obeyed to better than 0.5% in absolute terms, and to 
well with in the estimated errors in the exponents. Equa- 
tion (131) is obeyed to within 5%, a level consistent with 
the greater uncertainties in the values of f3 and 7. This 
agreement supports the validity of the scaling hypothesis 
presented in Eq. ([Tl]). 

Small-r expansion: There is an alternative method for 
determining the exponent x associated with the local spin 
susceptibility at the critical Kondo-coupling. The local 
spin-spin correlation function can be calculated pertur- 
batively in terms of po Jc ~ r i in analogy to the stan- 
dard e— expansion. Exponentiating the leading logarith- 
mic term in the perturbation series leads to 



1 - Jj(PoJc) 2 , 



PoJc < I- 



(16) 



Figure |l| compares this result with the NRG data for 
N = 2. The agreement is remarkably good, even when r 
(and hence poJ c ) is not small. 
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FIG. 1. Exponent 1 - x plotted vs r for TV = 2. The 
symbols plot the A = 9 data from Table [ij (Estimated er- 
rors are smaller than the symbolsj The solid line shows the 
perturbative prediction of Eq. Qlq) , in the form of a spline 
fit through numerical values for (poJ c ) 2 , calculated for A = 3 
and extrapolated to the continuum limit A = 1 (see Ref. jlq] ). 

For large but finite N, inserting the slave-boson mean- 
field value poJ c = r into Eq. ( |l6| ) yields x = 1 — 
2r 2 /N + 0(l/N) 2 . For N = oo, this expression repro- 
duces that obtained by combining slave-boson mean-field 
theory with the scaling ansatz [see Eq. (|i~5|)] . The agree- 
ment between the two results provides additional support 
for the scaling hypothesis. 

Finally, we note that the many-body spectrum at the 
critical point displays a non-Fermi liquid form. This is 
to be expected: since J c ^ 0, interactions do not renor- 
malize to zero. We have also studied the particle-hole- 
asymmetric problem, and again found the exponent x to 
be r-dependent. Details will be given elsewhere. 

In summary, we have analyzed the critical properties of 
the Kondo problem with a pseudogap in the conduction- 
electron density of states. We find that the temperature 
exponent of the local spin susceptibility at the critical 
point depends continuously on the power of the pseudo- 
gap. (For large N, such a dependence appears only at 
order 1/N.) The implications of our results for lattice 
heavy fermion systems can be seen within the dynam- 
ical mean- field approach j^o), which analyzes a lattice 
problem in terms of an impurity coupled self-consistently 
to a fermionic bath. The development of a pseudogap 
in the lattice system will lead to a pseudogap in the 
corresponding effective impurity problem. Should the 
self-consistency place the Kondo coupling near its crit- 
ical value, our results then imply that the local dynami- 
cal spin susceptibility will have an anomalous exponent. 
Tunneling and photoemission experiments in the relevant 
temperature regime would provide a direct test for the 
existence of a pseudogap. 
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